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Abstract According to various models, the orbital and 
the epicyclic frequencies of particles moving on a cir¬ 
cular orbit around compact objects are related to the 
quasi-periodic oscillations observed in the X-ray flux of 
some pulsars or black hole candidates. It is expected 
that they originate from the inner edge of the accretion 
discs, deep into the gravitational field of the compact 
objects. Considering the planned new generation X- 
ray timing observatories with large collective areas, the 
quasi-periodic oscillations might be an excellent tool for 
testing gravity in strong field regime and respectively 
alternative gravitational theories. We examine the or¬ 
bital and the epicyclic frequencies of a particle moving 
on a circular orbit around neutron or strange stars in B? 
gravity. The case of slow rotation is considered too. The 
gravity results are compared to the General Rela¬ 
tivistic case. We comment the deviations from General 
Relativity, as well as the deviations due to rotation in 
both theories. 

Keywords neutron stars ■ alternative gravity • 
oscillations 


1 Introduction 

In the last few years the interest in alternative theories 
of gravity was significantly increased. Major role for 
this has the experimental confirmation of the acceler¬ 
ated expansion of the universe and the fact that it does 
not ht in the predictions of General Relativity (GR) 
without the introduction of the so-called dark energy, a 
matter with exotic properties which interacts with the 
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visible matter only gravitationally. A class of viable al¬ 
ternative theories of gravity are the so-called f{R) the¬ 
ories. In this case we exchange the Lagrangian of the 
General Relativistic Einstein-Hilbet action, namely the 
Ricci scalar, with a more general one - a function of the 
Ricci scalar, hence f{R) theories . These theories 
allow us to resolve the problem with the accelerating 
universe without the necessity of introducing such an 
exotic matter. However, the theory should be tested 
not only on cosmological but on astrophysical scales 
too. The predictions, concerning astrophysical effects, 
should be very close to the GR ones in the weak field 
regime, but it is expected to deviate for strong fields. 
As a source of strong gravitational fields neutron stars 
(NS) and black holes (BH) are natural laboratories for 
testing alternative theories of gravity. In this work we 
are concentrating our efforts on the former. 

Future X-ray timing observatories with large col¬ 
lecting area like SKA Q, NIGER LOFT Q, and 
AXTAR @ may give us a chance to test gravity in 
strong field regime. A promising way to do that around 
compact objects are the so-called quasi-periodic oscil¬ 
lations (QPOs). QPOs are Hz to kHz oscillations in the 
X-ray flux of pulsars or black hole candidates. The lat¬ 
ter QPOs are supposed to occur only in the presence of 
strong gravitational fields, with their origin in the inner 
edge of the accretion disc, deep into the gravitational 
field of the star. However, the exact source of these os¬ 
cillations is unknown so there are different models for 
their explanation. In the beat frequency models, some 
connection between the orbital frequency and the spin 
frequency of the central compact object is suggested. 
These models require some azimuthally non-uniform 
structure co-rotating with the central object. In the rel¬ 
ativistic resonance models a resonance may occur at a 
particular radii of the disc at which the orbital and the 
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epicyclic frequencies have integer ratios. It is suggested 
that some form of resonance may occur due to the inter¬ 
action between the accretion disc and the central body 
too. The preferred radii models suggest that mecha¬ 
nism, that chooses some radii, exists. The relativistic 
precession models, investigated in some resent papers 
|9l4l3|. are based on the assumption that the QPOs are 
directly related to the orbital (Kepler) I7p and to the 
epicyclic frequencies. In the above-mentioned models 
all these frequencies or some of them are involved. The 
radial Vr and the vertical ve epicyclic frequencies occur 
if a particle on a stable circular orbit is perturbed. It 
will start to oscillate in a radial and in a vertical di¬ 
rection with some stable frequencies. This frequencies 
are the radial and the vertical epicyclic ones. A review 
of these models can be found in 1^. In some papers 
QPOs models concerning oscillations of the accretion 
disk itself are examined too iB- 17|. 


The epicyclic frequencies of rapidly rotating strange 
stars are examined in the resent paper [l^. Some of the 
above mentioned models were examined in alternative 
theories of gravity through the years [H, Ilfll - Elj l . The 
observed deviations from GR in these papers were quite 
low, except the case of rapidly rotating neutron stars 


in scalar-tensor theories [21 1. 


The structure of this paper is as follow: In Section 
II we present the basic steps for deriving the radius 
of the innermost stable circular orbit (ISCO), the or¬ 
bital frequency and the radial and the vertical epicyclic 
frequencies. In Section III we present and discus the 
results for neutron and strange stars. The paper ends 
with conclusions. 


2 The epicyclic frequencies of a rotating 
neutron star 


In this section we briefly present the basic steps in the 
derivation of the equations for the radius of the inner¬ 
most stable circular orbit (ISCO), the equations for the 
radial and for the vertical e picy clic frequencies and for 
the orbital frequency [l3) The equations de¬ 

scribing stable stationary neutron star models and ad¬ 
ditional mathematical details concerning gravity can 
be found in 25-27|. 

We are considering a stationary and axisymmetric 
spacetime with a metric 


ds^ = gudt^ + grrdr^ + geedO"^ -I- 2gt^dtdtp + g^^dip%l) 

where all the metric functions depend only on the co¬ 
ordinates r and 9. The massive particles subject to the 
gravitational force only move on timelike geodesics of 


the metric O- The stationary and axial Killing sym¬ 
metries of metric, generated by the Killing vectors ^ 
and give rise to two constants of motion, namely 
E = —ut and L = u^p. The first one corresponds to 
the energy per unit mass and the second one to the 
angular momentum along the axis of symmetry, and 
= dx^/dr is the four-velocity of the particle. 
It is not difficult one to show that the two conservation 
laws can be casted in the form 

_ ^g^p^p T ^9tip /rt\ 

d^ ~ ^2 ’ ^ ’ 

dip _ Egtp, + Lgtt . , 

dr g2 ' ^ ^ 

where we defined g2 = g^^ — gttg-pp for simplicity. Here 
t denotes the coordinate time, and r the proper time. 
From the normalization condition = —1, we 

also have 


grr'r'^ + ggg9^ + E‘^U{r,0) =-1 (4) 

with 


C/(r, 9 ) = 

52 

and I = LjE being the proper angular momentum.The 
derivatives in the above equations are with respect to 
the proper time r. 

For g = § the problem reduces to an effective one 
dimensional problem 

f^ = V{r), (6) 

with an effective potential 


y{r) = g^^ 


-l-E^U{r,9=^) 


( 7 ) 


The stable circular orbit with a radius f is deter¬ 
mined by the conditions V(re) = 0 = V (re) and V"{rc) > 
0, where with prime we denote the derivative with re¬ 
spect to r . The condition V"{rc) = 0 gives the ISCO 
radius. The angular velocity fip of a particle moving 
on a circular equatorial orbit can be found from the 
geodesic equations in the following way. We write down 
the geodesic equations in the form 


d f dx'^ \ I dx’^ dx'^ 
which for the radial coordinate gives 


( 8 ) 




2 a 


rgtp) 


dt dip ^ ^ ^ fdip 2 


dr dr 


drgepepi — 0 . 


( 9 ) 
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Taking into account that the angular velocity is defined 
by I7j, = ^ ^ we obtain from the above equation 

^ dip drgtip i \/ ipr^tip)^ drgttdrgipip /'1n^ 

To derive the epicyclic frequencies we should inves¬ 
tigate small perturbations of a stable circular orbit. The 
perturbations are written in the form 


r{t)=r + 5r{t), e{t) = ^ (11) 

where 5r{t) and 59{t) are perturbations to the stable 
circular orbit with coordinate radius f in the equatorial 
plane. The perturbations could be written explicitly in 
the form Sr{t) ~ S9{t) ~ Substitut¬ 

ing (na in eq. 0 and after some calculations, and a 
change of the proper time r with the coordinate one t, 
we obtain the expressions for the radial and the vertical 
epicyclic frequencies: 


{gtt + ^pgt^) 

2(27r) grr V 2/ 

(12) 

(Stt + c TTN ^ 

2(27r) gee V 2/ 

(13) 


For the case of static neutron stars the orbital fre¬ 
quency and the vertical epicyclic frequency coincide, i.e. 
V 9 = Vp, for / = 0, where / is the rotational frequency 
of the star (/ = ^). At the ISCO the square of the ra¬ 
dial epicyclic frequency is equal to zero, and for smaller 
radius it is negative, which shows a radial instabilities 
for orbits with radius smaller than the ISCO. 

For simple accretion disc models the inner edge of 
the disc is defined by the ISCO. For models with smaller 
masses the ISCO is in the interior of the star and for 
massive ones the radius of the ISCO is bigger than the 
radius of the star itself The inner edge of the accretion 
disc, therefore, can reach down to the surface of the star 
in the former case and to the ISCO in the latter one. 


3 Numerical results 


We investigate what are the changes in the radius of 
the ISCO and in the orbital and epicyclic frequencies in 
f{R) gravity with Lagrangian f{R) — the so- 

called R? gravity. The results are compared to pure GR. 
Deviations due to rotation, in slow rotation approxima¬ 
tion, are examined too. We consider two hadronic equa¬ 
tions of state (EOS) and a quark one. For the hadronic 
EOS we are using peacewise polytropic approximation 
281. The ones we choose are APR4 and MSI. APR4 


has maximal mass not much bigger than the observa¬ 
tional limit of two solar masses (medium stiffness), and 
MSI is a stiffer one with higher masses and radii. We 
examined a soft EOS too but the results are qualita¬ 
tively the same, and because of that, graphs are not 
presented. The quark EOS has the analytical form 


p = b{p-po), (14) 

where the constants b and po are taken from for 
EOS SQS B60. This EOS leads to maximal masses 
slightly bellow two solar masses, but we find it to be 
a good representative. 

In this section we examine static, f — 0 Hz, and 
slowly rotating, / = 80 Hz and / = 160 Hz, models 
of neutron and strange stars. Eirst we will discus the 
deviations from pure GR due to R^ gravity and as a 
next step we will examine the differences of the results 
for different rotational frequencies. The presented re¬ 
sults in all figures are up to the maximal mass. For the 
presented results in R^ gravity, the maximal deviation 
from GR is for the maximal adopted value of the pa¬ 
rameter a = 10^. This number is close to the maximal 
values allowed by the observations a ~ 10® or in dimen¬ 
sional units - a < 5 X [s^. 

In Fig. □ sequences of models, representing the ra¬ 
dius of the ISCO as a function of the stellar mass are 
plotted. A wide range of values for the parameter a is 
examined. In the left panel the results for EOS APR4 
are plotted, in the middle one for EOS MSI, and in 
the right one for EOS SQS B60. If the ISCO is in the 
interior of the star, that happens for models with low 
masses, we are plotting the radius of the star instead. 
The points in the different sequences of models, where 
the radius of the ISCO gets equal to the radius of the 
star are marked with an asterisk. For all examined EOS 
the results are qualitatively the same. Eor models with 
maximal masses the largest deviation is for the maxi¬ 
mal values of the parameter a. It is around 10 % and 
decreases with the decrease of a. In the limiting case of 
a —>■ 0 the solutions converge to the GR ones. 

At this point a comment concerning the presented 
results should be done. In Fig. [U in the present paper, 
the following behaviour for the models for which the 
radius of the ISCO gets bigger than the radius of the 
star can be observed. Let us start the discussion with 
the models marked with an asterisk. For small values 
of the parameter a they have smaller radius and lower 
mass than the GR case. The a = 10 models have param¬ 
eters which are close to the GR ones and for high values 
of a the radii and the masses are higher compared to 
the GR case. When the ISCO is outside the star its ra¬ 
dius is bigger compared to the GR case for models with 
equal masses and small values of a. For a = 10 the re- 
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suits are comparable to the GR ones. For bigger values 
of a the radius of the ISCO is smaller than the GR one. 
The explanation of this behaviour we find in the non- 
monotonous behaviour of some stellar parameters as a 
function of a. This is demonstrated, for example, in Fig. 
3 in 23 for the maximal mass. It was shown that for 


small values of a the mass decreases compared the GR 
one. At some value of the parameter a there is a turning 
point and the mass starts to increase, reaching maxi¬ 
mal masses higher than the GR one. Our investigations 
show that this is the case for neutron stars as well as 
for strange stars. Such behavior can be observed also 
for other stellar parameters, although in the different 
cases the minimum will be shifted to some other values 
of the parameter a. If one, for example, calculates the 
mass and the radius for models with fixed central en¬ 
ergy density in GR and in gravity, the same pattern 
can be found. It is not hard to see that there is a similar 
dependence for the models marked with an asterisk as 
well as for the case when the ISGO is outside the star. 
This leads us to the conclusion that this behaviour is 
characteristic one for the theory. Naturally, similar pat¬ 
terns can be found in the other graphs presented below. 


In Fig.[2]we plot the orbital frequency I'p = Qpl2TT, 
in kHz, at the ISGO as a function of the stellar mass. 
If the ISCO is in the interior of the star, we calculated 
Vp at the surface of the star. For all EOS we see qual¬ 
itatively the same results. The deviations from GR for 
models with a = 10^ is in average around 15 % and it 
does not change significantly with the mass. For mod¬ 
els with small value of a the deviation from GR, when 
the radius of the ISCO get bigger than the radius of 
the star, is close to the maximal one, but it rapidly de¬ 
creases with the increase of the mass, and converge to 
GR for maximal masses. 


In Fig. [ 3 ] we plot the maximal value of the radial 
epicyclic frequency, in kHz, as a function of the stellar 
mass. If the ISCO is in the interior of the star, we take 
the maximal value on the surface or outside the star. 
In the examined interval of masses the frequency mono- 
tonically decreases. For higher values of the parameter 
a the maximal radial epicyclic frequency is higher, com¬ 
pared to the GR one. For small values of a the frequen¬ 
cies are lower than the GR ones. For smaller masses 
there are higher deviations, but for maximal masses the 
plots converge to GR. 

The deviations due to rotation are qualitatively and 
quantitatively the same for the hadronic and for the 
quark EOS. This, combined with the small magnitude 
of the deviations is the reason why we chose to include 
only the favoured by the observations hadronc EOS in 
the following graphs. 


In Fig. H]we plot the difference between the orbital 
frequency i^p and the vertical epicyclic frequency vg, 
i.e. the nodal precession frequency t'„, as a function of 
the mass. In the left panel we plot the results for EOS 
APR4 and in the right panel for EOS MSI. The plots 
are for three different rotational rates: / = 0 Hz (the 
static case), / = 80 Hz, and / = 160 Hz. In black con¬ 
tinuous lines are the result for GR and in colour dashed 
line - the results for a = 0.3 and a = 10“^. As we said 
before, Vp and vg coincide in the static case. With the 
increase of the rotational frequency of the star the dif¬ 
ferences between the two frequencies increase.The dif¬ 
ference is calculated on the surface of the star if the 
ISCO is inside the star, and at the ISCO in the oppo¬ 
site case. The transition models are marked with as¬ 
terisks, in consistence with the previous graphs. The 
nodal precession frequency tends to decrease with the 
increase of the mass for models having ISCO in the in¬ 
terior of the star. It increases with the increase of the 
mass for models having ISCO outside the star. In the 
case of gravity we plot the results for a = 10'^ and 
a = 0.3. The former value of the parameter gives the 
maximal deviation from GR, and the latter we choose 
because of the interesting behaviour demonstrated in 
the graphs above. If the ISCO is in the interior of the 
star, the results are the same for all rotational rates. In 
the opposite case, when the ISCO is outside the star, 
the behaviour is the following. For a = 0.3, is slightly 
higher than the GR case, and rapidly converge to GR 
with the increase of the mass. For a = 10“^, has lower 
values than the GR one. The deviation from GR is more 
or less constant with the increase of the mass, and it is 
about 4 %. 


In Figs. O El [3 we plot the radius of the ISCO, the 
orbital frequency and the maximal radial epicyclic fre¬ 
quency as a function of the mass in the case of GR and 
the maximal deviation in R^ gravity (a = 10^) for the 
aforementioned values of /. The rotation have the same 
effect for all these quantities in GR and in R^ gravity. 
The reason why we are skipping the interestingly be¬ 
having case of a = 0.3 is the following. The changes 
of the ISCO and the examined frequencies for different 
rotational rates causes the graphs for a = 0.3 and a = 1 
to overlap with the GR ones. We find that case to be 
quite uninformative so only the maximal deviation is 
presented, namely a = 10^. Due to rotation the radius 
of the ISCO decreases with about 4 %, and the frequen¬ 
cies increase with the same magnitude. The percentage 
deviations are the same for GR and for R^ gravity. 
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Fig. 1 The radius of the ISCO in km as a function of the stellar mass for different values of the parameter a. In the left panel 
the results for EOS APR4 are plotted, in the middle one for EOS MSI, and in the right one for SQS B60. The results are for 
static neutron stars (f = 0 Hz). 



Fig. 2 The orbital frequency of a particle on a circular orbit in kHz as a function of the mass of the star for different values 
of the parameter a. In the left panel the results for EOS APR4 are plotted, in the middle one for EOS MSI, and in the right 
one for SQS B60. The results are for static neutron stars (f = 0 Hz). 


4 Conclusions 

In this paper we examined the orbital and the epicyclic 
frequencies of particles moving in circular orbits around 
neutron and strange stars. Various models relate all 
these frequencies or some of them to QPOs. Taking 
into consideration the new generation observatories for 
timing of X-ray pulsars with large collective area which 
are under construction now, QPOs are expected to have 
major role in testing strong field regime of gravity. 

gravity is a viable alternative to GR. We exam¬ 
ined how changing the parameter of the theory effects 
the radius of the ISCO, the orbital and the epicyclic fre¬ 
quencies. Most of the observed neutron stars are slowly 
rotating, so we concentrated our efforts on the static 
and slowly rotating solutions. We investigated the de¬ 


viations of the slow rotating solutions from the static 
case too. 

We concerned two hadronic and one quark equations 
of state and for both cases we observed qualitatively the 
same behaviour with the change of the parameter of the 
theory. The largest deviations from GR occurs for the 
maximal adopted value of the parameter a. The radius 
of the ISCO decreases with about 10% and the orbital 
and the maximal radial epicyclic frequencies increases 
with 15 - 20 %. However, the effect due to slow rotation 
is expectedly small. For the highest examined rotational 
frequency (/ = 160 Hz) the deviation from the static 
case is about 4% for all examined quantities. 

The deviations for slowly rotating stars, up to 160 
Hz, are much smaller than the deviations due to the 
gravity. Because of that static theoretical and numerical 
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Fig. 3 The maximal value of the radial epicyclic frequency in kHz as a function of the stellar mass for different values of the 
parameter a. In the left panel the results for EOS APR4 are plotted, in the middle one for EOS MSI, and in the right one for 
SQS B60. The results are for static neutron stars (f = 0 Hz). 



M/M0 M/M0 


Fig. 4 The difference between the orbital and the vertical epicyclic frequency at the ISCO. The plots are for different rotational 
rates for GR and for gravity with a = 0.3 and a = 10^. In the left panel the results for EOS APR4 are plotted, and in the 
right one for EOS MSI. 



Fig. 5 The radius of the ISCO, in km, as a function of the stellar mass. The plots are for GR and for gravity with a = 10'^ 
and for stars rotating with / = 0 Hz, / = 80 Hz, and / = 160 Hz. In the left panel the results for EOS APR4 are plotted, and 
in the right one for EOS MSI. 
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Fig. 6 The orbital frequency at the ISCO, in kHz, as function of the stellar mass. The plots are for GR and for B? gravity 
with a = 10^ for stars rotating with / = 0 Hz, / = 80 Hz, and / = 160 Hz. In the left panel the results for EOS APR4 are 
plotted, and in the right one for EOS MSI. 



Fig. 7 The maximal radial epicyclic frequency as a function of the stellar mass. The plots are for GR and for gravity with 
a = 10^ for stars rotating with / = 0 Hz, / = 80 Hz, and f = 160 Hz. In the left panel the results for EOS APR4 are plotted, 
and in the right one for EOS MSI. 


models can be used for comparison with observational 
data even for slowly rotating stars. 

The frequencies examined in this paper provide us 
with an opportunity to test GR in strong field regime, 
as well as alternative theories. Results from the numeri¬ 
cal models could be compared to the observational data 
expected from the new generation X-Ray timing obser¬ 
vatories. However, even if the observational data do not 
provide a concrete information about the correct grav¬ 
itational theory, it will provide us with an unique op¬ 
portunity to restrict the wide range of possible values 
for the free parameter of the theory. 
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